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Abstract 



It is well-known that one-dimensional isentropic gas dynamics has two elementary 
waves, i.e., shock wave and rarefaction wave. Among the two waves, only the rar- 
efaction wave can be connected to vacuum. Given a rarefaction wave with one-side 
vacuum state to the compressible Euler equations, we can construct a sequence of 
solutions to one-dimensional compressible isentropic Navier-Stokes equations which 
converge to the above rarefaction wave with vacuum as the viscosity tends to zero. 
Moreover, the uniform convergence rate is obtained. The proof consists of a scaling 
argument and elementary energy analysis based on the underlying rarefaction wave 
structures. 

Keywords: compressible Navier-Stokes equations, zero dissipation limit, rarefaction 
wave, vacuum. 



1 Introduction and main result 

In this paper, we investigate the zero dissipation limit of the one- dimensional com- 
pressible isentropic Navier-Stokes equations 




= e u 
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where p(t, x) > 0, u(t, x) and p represent the density, the velocity and the pressure of the 
gas, respectively and e > is the viscosity coefficient. Here we assume that the viscosity 
coefficient e is a positive constant and the pressure p is given by the 7— law: 

P(p) = P ~ 

with 7 > f being the gas constant. 

Formally, as e tends to zero, the limit system of the compressible Navier-Stokes equa- 
tions fll.ip is the following inviscid Euler equations 

(h + (pu)x = 0, 

(1-2) 

(pu) t + (pu +p(p)) x = 0. 

The Euler system (jl.2p is a strictly hyperbolic one for p > whose characteristic fields 
are both genuinely nonlinear, that is, in the equivalent system 



P + U P ) P 

u I \ p'(p)/p u I \ u 



0. 



t 



the Jacobi matrix 

u p 



p'(p)/p u 



has two distinct eigenvalues 



Ai(p,u) = u- y/tfip), A 2 (p,tt) = u + ^p'(p) 
with corresponding right eigenvectors 

r i (p ) «) = (l,(-l)*^S)* j , = i,2, 

P 

such that 

nip, u) ■ V P , U \( P , u) = ( _ ir ^>) + V( p) ^ Qj . = X) 2 

We can define the i— Riemann invariant (i = 1, 2) by 



S J (p,n) = n + (-ir +1 | P ^S3 

such that 



s 



V(p,u)^i(p,u) ■ n(p, u) = 0, Vp>0,w. 
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The study of the limiting process of viscous flows when the viscosity tends to zero, is 
one of the important problems in the theory of the compressible fluid. When the solution 
of the inviscid flow is smooth, the zero dissipation limit can be solved by classical scaling 
method. However, the inviscid compressible flow contains singularities such as shock and 
the vacuum in general. Therefore, how to justify the zero dissipation limit to the Euler 
equations with basic wave patterns and/or the vacuum is a natural and difficult problem. 

There have been many results on the zero dissipation limit of the compressible fluid 
with basic wave patterns without vacuum. For the system of the hyperbolic conservation 
laws with artificial viscosity 

u t + f(u) x = eu xx , 

Goodman- Xin [5] first verified the viscous limit for piecewise smooth solutions separated 
by non-interacting shock waves using a matched asymptotic expansion method. Later 
Yu [22] proved it for the corresponding hyperbolic conservation laws with both shock 
and initial layers. In 2005, important progress made by Bianchini-Bressan[T] justifies the 
vanishing viscosity limit in BV space even though the problem is still unsolved for the 
physical system such as the compressible Navier-Stokes equations. For the compressible 
isentropic Navier-Stokes equations ( 11 .11) . Hoff-Liu [6] first proved the vanishing viscosity 
limit for piecewise constant shock even with initial layer. Later Xin [20] obtained the zero 
dissipation limit for rarefaction waves without vacuum for both rarefaction wave data and 
well-prepared smooth data. Then Wang [IB] generalized the result of Goodmann-Xin [5] 
to the isentropic Navier-Stokes equations ( 11. ip . For the full Navier-Stokes equations where 
the conservation of the energy is also involved, there are also many results on the zero 
dissipation limit to the corresponding full Euler system with basic wave patterns without 
vacuum. We refer to Jiang-Ni-Sun [11] and Xin-Zeng [21] for the rarefaction wave, Wang 
[19] for the shock wave, Ma [14] for the contact discontinuity and Huang- Wang- Yang 
[9j [10] for the superposition of two rarefaction waves and a contact discontinuity and the 
superposition of one shock and one rarefaction wave cases. 

More recently, Chen-Perepelitsa [3] proved the vanishing viscosity to the compressible 
Euler equations for the compressible Navier-Stokes equations (II. ip by compensated com- 
pactness method for the general case if the far field of the initial values of Euler system 
( II. 2p has no vacuums. Note that this result is quite universal since the initial values of the 
Euler system can contain vacuum states in the interior domain. Huang-Pan- Wang- Wang- 
Zhai [7] established the corresponding results to the compressible Navier-Stokes equations 
( II. ip with density-dependent viscosity. 
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Now we turn back to the case of the basic wave patterns with vacuum states. As pointed 
out by Liu-Smoller [13], among the two nonlinear waves, i.e., shock and rarefaction waves, 
to the one-dimensional compressible isentropic Euler equations (II .21) . only the rarefaction 
wave can be connected to vacuum. However, to our knowledge, so far there is no any 
results on the zero dissipation limit of the system (II .ip in the case when the Euler system 
(II. 2p contain the rarefaction wave connected to the vacuum. In this paper, we investigate 
this fundamental problem and want to obtain the decay rate with respect to the viscosity 
e. Remark that Perepelitsa [T6j consider the time-asymptotic stability of solutions of 1-d 
compressible Navier-Stokes equations ( 11.11) toward rarefaction waves connected to vacuum 
in Lagrangian coordinate and Jiu-Wang-Xin [12] study the large time asymptotic behavior 
toward rarefaction waves for solutions to the 1-dimensional compressible Navier-Stokes 
equations (II. ip with density- dependent viscosity for general initial data whose far fields 
are connected by a rarefaction wave to the corresponding Euler equations with one end 
state being vacuum. 

Now we give a description of the rarefaction wave connected to the vacuum to the 
compressible Euler equations ( 11.21) . see also the references [13] and [17]. For definiteness, 
2-rarefaction wave will be considered. If we investigate the compressible Euler system 
(II. 2p with the Riemann initial data 



where the left side is the vacuum state and p + > 0, u + are prescribed constants on the 
right state, then the Riemann problem (11.21) . (11.31) admits a 2— rarefaction wave connected 
to the vacuum on the left side. By the fact that along the 2— rarefaction wave curve, 
2— Riemann invariant £ 2 (p, u) is constant in (x, t), we can get the velocity w_ = £ 2 (p + , u + ) 
being the speed of the fluid coming into the vacuum from the 2-rarefaction wave. This 
2— rarefaction wave connecting the vacuum p = to (p + ,w + ) is the self-similar solution 



p(0,x) =0, x < 0, 
(p,u)(0,x) = (j>+,u+), x > 0, 



(1.3) 



(p 



,u r2 ){£), (f = f) of defined by 




p r2 (O = 0, if £<A 2 (0, 
f, if u- < £ < A 2 (p+,u + ), 

A 2 (p+,w+), if £ > A 2 (p+,u+), 



(1.4) 



and 



wm« r2 (0) = s 2 (o,«_) = s 2 (p + , M+ ). 



(1.5) 
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Thus we can define the momentum of 2-rarefaction wave by 




0, if p r2 = 0. 



(1.6) 



In the present paper, we want to construct a sequence of solutions (p e ,m £ )(x,t) to the 
compressible Navier-Stokes equations ( 11.11) which converge to the 2-rarefaction wave 
(P2, m r 2 ){x/t) defined above as e tends to zero. The effects of initial layers will be ignored 
by choosing the well-prepared initial data depending on the viscosity for the Navier-Stokes 
equations. 

The main novelty and difficulty of the paper is how to control the degeneracies caused 
by the vacuum in the rarefaction wave. To overcome this difficulty, we first cut off the 
2-rarefaction wave with vacuum along the rarefaction wave curve. More precisely, for any 
p > to be determined, the cut-off rarefaction wave will connect the state (p, it) = (p, u^) 
and (p+,zt+) where can be obtained uniquely by the definition of the 2-rarefaction 
wave curve. Then an approximate rarefaction wave to this cut-off rarefaction wave will be 
constructed through the Burgers equation. Finally, the desired solution sequences to the 
compressible Navier-Stokes equations fll.il) could be established around the approximate 
rarefaction wave. The uniform estimates to the perturbation of the solution sequences 
around the approximate rarefaction wave can be got by the following two observations. 
One is the fact that the viscosity e can control the degeneracies caused by the vacuum in 
rarefaction waves by choosing suitably p = p(e). In fact, we choose p = e a \ lne| with a 
defined in (13.111) in the present paper. The other observation is that we can carry out the 
energy estimates under the a priori assumption that the perturbation is suitably small 
in H 1 (R,) norm with some decay rate with respect to e as e tends to zero. See (13.101) in 
the below for the details. Note that this a priori assumption is natural but is first used 
in studying zero dissipation limit to our knowledge. With these two observations, we can 
close the a priori assumption and obtain the desired results. 

Now our main result is stated as follows. 

Theorem 1.1. Let (p r2 , m T2 )(x/t) be the 2-rarefaction wave defined by (ll.4jl - fll.6p with 
one-side vacuum state. Then there exists a small positive constant €q such that for any 
e G (0,eo), we can construct a global smooth solution (p e ,m e = p € u e )(x,t) with initial 
values (13. ip to the compressible Navier-Stokes equation U.l\) satisfying 



(V 



(p e - p T \nf - m r2 ), (p% m% E C°((0, +oo); L 2 (R)) 

rrf xx e L 2 (0,+oo;L 2 (R)). 
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2) As viscosity e — > 0, (p e ,m e )(x,t) converges to (p T2 ,m T2 )(x/t) pointwisely except the 
original point (0, 0). Furthermore, for any given positive constant h, there exists a constant 
Ch > 0, independent of e, such that 

sup||p e (-,t)-p r2 (-)|| LOO <C ft e a |lne|, 

t>h t 

II 't A r 2( -M\ ^ J Ch*\lne\-t, ^ 1 < 7 < 3, ( L7 ) 
sup \\m (•, t) — m ^(-) kao < < 1 

t>h V ' V t 7 " ~ \ C h e— |lne|, if 7 > 3, 

with the positive constants a, b given by 

if 1 < 7 < 2, 

(1.8) 

if 7 > 2. 



and 




4(7 + 4) 



if 1 < 7 < 2, 

(1.9) 

if 2< 7 <3. 



A few remarks are followed. 



Remark 1.2. Similar result to Theorem \l.l\ is also expected for a one- dimensional com- 
pressible Navier-Stokes equation with density- dependent viscosity which reads 

Pt + (pu)x = 0, 

(pu) t +(pu 2 + p^ x = e(p a u x ) x , (1 ' 10) 

with suitable a > and 7 > 1. Actually, the system (11.101) can be derived by Chapman- 
Enskog expansions from the Boltzmann equation where the viscosity of the compressible 
Navier-Stokes equations depends on the temperature and thus on the density for isentropic 
flows. Also, the viscous Saint-Venant system for the shallow water, derived from the in- 
compressible Navier-Stokes equation with a moving free surface, is expressed exactly as 
in (ll.lOp with a = 1 and 7 = 2, see [H [^]. In this situation, since viscosity vanishes at 
vacuum, the convergence rate with respect to e may become slower than in Theorem \1.1\ 
and may depend on a and 7. However, this is left to the forthcoming paper. 

Remark 1.3. Our result and method can also be generalized to the 1-D full compressible 
Navier-Stokes equations with the conservation of the energy and the Boltzmann equation 
with slab symmetry. This is left to the forthcoming paper. 
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Remark 1.4. It is also interesting to study the zero dissipation limit of compressible 
Navier-Stokes equations U.l\) in the case when the Euler system U.ty) has two rarefaction 
waves with the vacuum states in the middle. However, it is nontrivial to cut off these 
rarefaction waves with vacuum along the corresponding rarefaction wave curves. In fact, 
the wave structure containing two rarefaction waves with the medium vacuum is destroyed 
and some new wave may occur in the cut-off precess, which is quite different from the 
single rarefaction wave case considered in the present paper. 

Remark 1.5. It is noted that in the a priori estimates {3.11) below , the estimates for cf) 2 
from the potential energy holds with the weight p 7 ~ 2 which is degenerate at vacuum when 
7 > 2. Therefore, the convergence rate obtained in Lemma WTB and thus in Theorem li.il 
depends on 7 when 7 > 2. 

The rest of the paper is organized as follows. In section 2, we construct a smooth 
2-rarefaction wave which approximates the cut-off rarefaction wave based on the inviscid 
Burgers equation. And the proof the Theorem ll.il is given in Section 3. 

Throughout this paper, H l (H), I = 0,1, 2, denotes the l-th order Sobolev space with 
its norm 

1 

ll/lh = (£lW)*' andll-HHHUw, 
3=0 

while L 2 (dz) means the L 2 integral over R with respect to the Lebesgue measure dz, 
and z = x or y. For simplicity, we also write C as generic positive constants which are 
independent of time t and viscosity e unless otherwise stated. 



2 Rarefaction waves 

Since there is no exact rarefaction wave profile for the Navier-Stokes equations (11.11) . 
the following approximate rarefaction wave profile satisfying the Euler equations was 
motivated by Matsumura-Nishihara [15] and Xin [20]. For the completeness of the presen- 
tation, we include its definition and the properties listed in Lemma 12.11 However, Lemma 
12. II is a little different from [20] as stated after Lemma [2.11 

Consider the Riemann problem for the inviscid Burgers equation: 

w t + ww x = 0, 

w_, x<0, (2.1) 



w(x, 0) 



w + , x > 0. 
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If w^ < w + , then the Riemann problem (12 .ip admits a rarefaction wave solution w r (x, t) 
w r (j) given by 



™ r (§)=< ?, w^<^<w +J (2.2) 



As in [20], the approximate rarefaction wave to the compressible Navier-Stokes equations 
(II. ip can be constructed by the solution of the Burgers equation 

w t + ww x = 0, 

x w + + w- w+-w- x (2.3) 
w{0,x) = ws{x) = w{-) = 1 tanh-, 

OA 2 

where 5 > is a small parameter to be determined. In fact, we choose 6 = e a in ( 13.121) 
with a given by (13.111) in the following. Note that the solution Wg(t,x) of the problem 
(12. 3p is given by 

Wg(t,x) =ws(x (t,x)), x = x Q (t,x) + ws(x Q (t,x))t. (2.4) 
And Wg(t,x) has the following properties: 

Lemma 2.1. The problem <\2.3§ has a unique smooth global solution Wg(x,t) for each 
5 > such that 

(1) w- < w r 5 (x,t) < w + , d x w r s {x,t) > 0, for x e R, t > 0, 5 > 0. 

(2) The following estimates hold for all t > 0, 5 > and pE [1, oo]: 

\\d x w r 5 (;t)\\ LP <C(w + -w-) 1 ' p (5 + t)- 1+1 '>, (2.5) 

\\d 2 x w r s (;t)\\ LP <C(5 + t)- l 5- 1 ^^ (2.6) 

^ 2 w r s (x,t) < A dw r s (x,t) 
dx 2 ~ 5 dx 

(3) There exist a constant <5 G (0, 1) such that for 5 G (0, So], t > 0, 

\\w r s (-,t) - w r (-)\\ L ~ < C5t~ l [ln(l + t) + | lntf|]. 



S 



The proof of Lemma [27T1 can be found in Xin [20]. However, the description of Lemma 
12.11 is equivalent to but a little different from Xin [20]. Take the estimation (I2.5P as an 
example, which is described by 

\\d x w r s (-,t)\\ LP < Cmin{K - w_)5- 1+1 /v, (w + - w^r 1+1 ^}, (Q' 

in Xin [20] • In fact, two estimations f !2.5p and (12. 5p ' are equivalent for fixed wave strength 
w + — w^. However, the advantage of Lemma 12.11 is that the energy estimate can be 
carried out for all time since there is no singularity to the approximate rarefaction wave 
even at t = 0. While in Xin's paper [20], the energy estimate must be done in two time- 
scalings, that is, finite time and large time, due to the singularity of the estimations of 
the approximate rarefaction wave at t = 0. 

As mentioned in the introduction, we will cut off the 2-rarefaction wave with vacuum 
along the wave curve in order to overcome the difficulty caused by the vacuum,. More 
precisely, for any p > to be determined, we can get a state (p, u) = (p, u At ) belonging to 
the 2-rarefaction wave curve. From the fact that 2-Riemann invariant £ 2 (p, u) is constant 
along the 2-rarefaction wave curve, can be computed explicitly by = £ 2 (p+,w+) + 
;pj-p 2 2 _ . Now we get a new 2-rarefaction wave (p 7 ^ ,u r ^)(^), (£ = x/t) connecting the 
state (n,Ufj) to the state (p + ,w + ) which can be expressed explicitly by 

{A 2 (p,?v), £<A 2 (p,?v), 
£, A 2 (p,u M )<£<A 2 (p + ,u + ), (2.8) 

A 2 (p+,^+), £>A 2 (p + ,u + ). 

and 

s 2 (p; 2 , M ; 2 ) = e 2 (p, Mm ) = s 2 ( P+ , M+ ). (2.9) 

Correspondingly, we can define the momentum function m^ 2 = p r ^u r ^. It is easy to show 
that the cut-off 2-rarefaction wave (p^ 2 , m r ^)(x/t) converges to the original 2-rarefaction 
wave with vacuum (p r2 , m T2 )(x/t) in sup-norm with the convergence rate p as p tends to 
zero. More precisely, we have 

Lemma 2.2. There exist a constant po £ (0, 1) such that for p G (0,po],t > 0, 

II (#,"»?)(•/*) - (p r2 ,m- 2 )(-/t)||^ < Cp. 
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The proof of Lemma 12.21 can be obtained directly from the explicit solution formula 
of rarefaction waves, so we omit it for brevity. 

Now the approximate rarefaction wave (p^s, u^s){x, t) of the cut-off 2-rarefaction wave 
)(f) t° compressible Navier-Stokes equations (II. ip can be defined by 

w+ = A 2 (p+, u+), = A 2 (p, u M ), 

w r s (t,x) = \ 2 (p^5,u^ s )(t,x), (2.10) 
^(p^SjU^sliXjt) = £ 2 (p + ,u + ) = S 2 (/i, Up), 

where w r s is the solution of Burger's equation (12. 3p defined in (I2.4|) . From then on, the 
subscription of (ps tfl ,us tfl )(x,t) will be omitted as (p,u)(x,t) for simplicity. Then the ap- 
proximate cut-off 2-rarefaction wave (p, u) defined above satisfies 

Pt + (pu) x = 0, 

/ o \ (2.H) 

(pu) t +(pu 2 +p(p)) x = 0, 

and the properties of the approximate rarefaction wave (p, u) is listed without proof in 
the following Lemma. 

Lemma 2.3. The approximate cut-off 2-rarefaction wave (p,u) defined in ( 12 .101) satisfies 
the following properties: 



(i) u x (x, t) = ^{w r s ) x > 0, for x G R, t > 0; 
Px = P^u x , and p xx = p^Uxx + ^p 2 ~ 7 (u 



(ii) The following estimates hold for all t > 0, 5 > and pG [1, oo]: 

\\u x (-,t)\\Lp<C{w + -w-) l /P(5 + t)- 1+1 /P, 

\\u xx (;t)\\ LP <C(5 + t)-M- 1+1 /P. 

(Hi) There exist a constant 5 G (0, 1) such that for 5 G (0, 5 ],t > 0, 

|| (p p^ 2 , u t) Hioo < C6t- 1 [ ln(l + t) + | In 5|] . 

3 Proof of Theorem [LI] 

To prove Theorem ll.il the solution (p e ,u e ) is constructed as the perturbation around 
the approximate rarefaction wave (p,u). Consider the Cauchy problem for (II. ip with 
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smooth initial data 

(p> e )(x,t = 0) = (p,n)(a;,0). (3.1) 
Then we introduce the perturbation 

(0, i[>)(y, r) = (p e , u%x, t) - (p, u){x, t), (3.2) 

where y, r are the scaled variables as 

2/ = ~, T = \' ( 3 - 3 ) 

and (p e , u € ) is assumed to be the solution to the problem (II. ip . For the simplicity of the 
notation, the superscription of (p e , u e ) will be omitted as (p, w) from now on if there is no 
confusion of the notation. Substituting ( 13.21) and ( 13.31) into ( II. ip and using the definition 
for (p,u), we obtain 



4>r + P?Py + U(f) y = -/, 
P^r + pWfjy + P'(p)(f)y ~ i)yy = 

(M)(y,0) = 0, 



where 



f = u y (j) + pyi/j, 

g = -Uyy + P^Uy + ^ 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



We seek a global (in time) solution (<f>, ip) to the reformulated problem ( 13.41) — ( 13. 6p . 
To this end, the solution space for ( 13. 4 p — ( 13 .6p is defined by 



with < T\ < +oo. 



» (^V)6G°([0 ) ti];^ 1 (R)) 



GL 2 (0,r i; L 2 (R)), 



Theorem 3.1. The problem flff.^P — (\3.6i) admits a unique global-in-time solution (0, ip) G 
X(0, +oo). Furthermore, there exist positive constants e and C independent of e, such 
that if < e < eo , i/ien 



sup 

rG[0,+oo] JR 

- + 0O 



j y l p 2 + pr- 2 <f ) 2 + <t> 2 y + ^)(T,y)dy 



+ 



R 



2 



^ + p 7 " 2 ^ 2 + pu^ 2 + (P-^I + dydr < Ce^ 2 ~^\ lner 1 



/2 



(3.? 
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Consequently, 



Ce^llnel" 1 / 4 , if 1< 7 < 2 

sup 0(-,t)l~<< 1 (1 V4 
o<r<+oo CW 4 |lne| l 7J/ , if 7 > 2, 



SUp ||^(-,T)||ioc < 



(3.9) 

Ce^llnel" 1 / 3 , if 1< 7 < 2, 



^vAjj 11 y v 5 ' / 1 1 — ■ : ^ 7+1 1 1 /o 

0<T< + OO Ce 4 ^+ 4 ) |lne|~ i/z , if 7 > 2. 

In what follows, the analysis is always carried out under the a priori assumptions 



sup U(;t)\\ Loo < e a , sup Uy\\<l, (3.10) 

0<t<ti r£[0,Ti] 



with a given by 



«> 1<7<2, 

6 1 (3.11) 

7 + 4 

Take 

/i = e a |lne|, 5 = e a , (3.12) 

in the sequel. Then it follows that \i > 2e a if e 1. Under the a priori assumption (I3.10p . 
we can get 

l<P<f (3-13) 

In fact, if e 1, then 

P = P + > P - > P - e a > p - \\i > (3.14) 

P = P + < P + < P + e a < P + \p < y • (3.15) 

Moreover, under the a priori assumption ( 13 .101) . it holds that 

Ci^V < p(p) - p(p) - p'(p)<f> < C 2 (P- 2 <j> 2 . (3.16) 

where Ci,C 2 are positive constants independent of e. 

Since the proof for the local existence of the solution to (13. 4p — (13. 6 p is standard, we 
omit it for brevity. To prove Theorem 13.11 it is sufficient to prove the following a priori 
estimates. 
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Lemma 3.2. (A priori estimates) Let 7 > 1 and (<p,ip) € X(0,Ti) be a solution to 
the problem ( 3.4\ ) — ( \3.6i) . Then under the a priori assumption (I3.10p . there exist positive 
constants eo and C independent of e, such that ifO<e<e , then 



re[0,-ri] Jr. 

fTl 



( pi) 2 + f~ 2 <p 2 + <f>l + il>l) (r, y)dy 



+ 



JR 



dydr^Ce^-^llne]- 1 / 2 . 



(3.17) 



Consequently, 



sup 

0<T<T1 



T Too < 



Ce 1/6 |lne|- 1/4 , if 1< 7 < 2, 

Ce^4|lne| (1 " 7)/4 , if 7 > 2, 

if 1 < 7 < 2, 



Ce^llne!" 1 / 2 , 
sup |M-,t)||£» < \ __2±i M ,_ 1/2 
o<r<n Ce 4 ^+ 4 ) I lne| v , if 7 > 2. 



(3.18) 



Proof of Lemma 13.2b The proof of Lemma 13.21 consists of the following steps. 
Step 1. First, define 

^ 2 



where 



E:= $(p,p) + 



Direct computations yield 



(3.19) 



pE) + puE - ipyip + {pip) - p{p))i) 

/ T L J y 

+*Py + «v (p(p) - 'P(P) - P'(P)<P) + = Uyyi). 

Then integrating the above equation over R 1 x [0,r] and using (13.131) . (13 . 1 61) and (13. 191) 
imply 

pip 2 + p 7 ~ 2 (f) 2 ^jdy + I I (^p 2 + p J ~ 2 u y (f) 2 + puyi/j^jdydr < C J J \u yy ip\dydT. 

(3.20) 



</R 
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By Sobolev inequality and Lemma I2.3[ one has 

\u yy i)\dydr<c I ||uJUi|Ml 1/2 ||^J 1/2 dT 

Jo 

1 




JR 



< 


c 




1 


< 


8 




1 


< 


8 




1 


< 





r + 5/e 

W^fdr + C I ( 



-) 4/3 ||^H 2/3 ^ 
r + d/e 

dr + - snp ll^f + Cf/i" 1 / 3 / (— t^) 4 ^ 

8 T6[o,n] V Jo T + o/e 

fm\ 2 dr+ l - snp ||^H 2 + C(^) 1/2 . 
Jo 8 r6[o,n] P d 

Combining ( I3T20D and (13T211 and recalling ( I3TT2D yield 



p^ 2 + f 2 2 ) (r, + J Jjfi + P 1 + P%^ 



(3.21) 



3/2 



sup 

re[0,n] ^R 



dydr 



(3.22) 



<C e (V2-a)| lne |-l/2. 

SYep We make use of the idea in [8] with modifications to derive the estimation 
of 4> y . Differentiating (13. 4p with respect to y and then multiplying the resulted equation 

by (f> y /p 3 to get 



2p 3 - 



v 2 P 3 >y 

Multiplying (13. 5p by <p y / p 2 gives 



l(Uyy(p + Pyy^ + 2p y ^y). 



P 



(3.23) 



-u n , 



p p p y p l p l 

.-. ^ , 9 - , * ^ 2 , ^ „- ^0 __ ^ 

T ^Py~ 



(3.24) 



P 



+ Pw— + PyUy-p ~ (My ^ 



p 2 



Adding (I3.23P and ( I3.24[) together, then integrating the resulted equation over R 1 x [0, r] 
imply 




<j>y H^y 
2p 3 P 

2 



dy 




V\p)-^dydr 



P 



'0 JR. 



^2/ + 2 P: 



P 



Pw P^ M 2/ p2 P M S/ ^2 



(3.25) 



__ f + Pj/y^ + 2pj,-0y) - g^^dydr. 
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The combination of (I3.22p and (13.251) leads to 



| + + J?- 2 <A dy+ [ [ U 2 y + p^uycf + pu v ij 2 + p^ci) 2 ) dydr 

R V P 7 JO JR V 



<C [ [ {K^I + \p y ^p-\ + \pyy%~\ + \pyUy^-\ + + \py^y 



! j n L ■ - p ■ ' p p- ■ V " /> V " 



^01 + 1%^. 
p A p A p 



\ g ^\\dydr + Ce^ 2 -^\]ne\-^ 2 



:= ^ /i + C < e (l/2- a )| lne | 



-1/2 



(3.26) 

Now we estimate the terms on the right hand side of (13.261) one by one. By Lemma 12.31 
(I3.13P and Cauchy's inequality, it holds that 



h = \u y ^\dydr<- U y \\ z dT + C / p^Uy^^-dydr 

Jo Jr P ° Jo Jo Jr P ' 

< I [ Uyfdr + Cp-^S^e [ [ p^- 2 u y( j) 2 dydT (3.27) 

'0 JR 



< \f Uvfdr + \ jf j^- 2 u y <p 2 dydr 



'0 ° JO JR 

where we have used the fact that 

Cp-^5~ l e = Ce^^llnel" 7 < Ce^|lne|- 7 < -, if e < 1. (3.28) 

8 

From Lemma [2.31 (i). one has 

p y = P^Uy, (3.29) 

one can get 

h = ^ I \Py^\dydT 

Jo Jr P 



— ~ I ll^y l| 2 ^ r + C / I P y u y {puyip 2 )dydr 




o Jr 



(3.30) 



<- I U y \\ 2 dT + Cp^6- l e [ [ puyi) 2 dydr 

H T i r f 

~8J ^y^ dT+ sj Q J^pUy^dydr 
where in the last inequality we used the fact (13. 28ft . 

Recalling (12 .7p from Lemma [2.11 and the fact (i) in Lemma [2 .3[ one can arrive at 

I- I I 2 _3^7 r 2(3 — 7)_o_- v/ r n2 i 

1^1 = \— x p * c fe + ^ i 3 

<c(^| + ^»f 
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Thus one has 
h 

< Ce5~ l 



nip 2 
\ Pyy —\dydT 
_L P 

nlpUy^p'^ldydr + CeS' 1 / / \pu y ^P''p' 1 \dyd'\ 
_t Jo Jr 



< C/i 7 5 1 e I I pu y tp 2 dydT 
Jo Jr 



< - I I pu y ip 2 dydr, if f < 1. 
8 Jo Jr 

It follows from Lemma 12.31 and (I3.29|) that 

\pyUy^-\dydr 
lo Jr. P 




- q I I PUyTp 2 dydr + C [ [ Uyf 2 <p 2 -^-dydr 
° Jo Jr Jo Jr. P 



-\f o f R pUy ^ 2dydT + c j^f§2 f f R pl 2u v ( t )2d y dT 

-If I puy^dydr + - [ I f^Uy^dydr. 
° Jo Jr. ° Jo Jr 



Similarly, it holds that 



h = I I \uy^—^-\dydr 
'o Jr P 



< I rWP^U'dr + C [ T [ pUy^dydT 
° Jo Jo Jr P 1 

~\So \\ p2 ^^ 2dT + C J^ J J^v^dydr 
-q I IIP^<M 2rfr + 1 I I pu y ip 2 dydr. 



° JO ° JO JR 

By Lemma the equality (I3.29P and Cauchy's inequality, one has 




\py^ y \dydr 
'o Jr P 



< \ f WP^U^t + C [ T [ %^ 2 ydydr 
° Jo Jo Jr P 



< 



11 ^ (7 ~ 3)/2 ^ 2 ^ + ^/ L^ dvdr 



<\ f\\p { ^ )/2 U 2 dr+\ I I ridydr. 



1(3 



y 

'0 JR 



Similarly, I7 can be estimated as 

, u 



h 




Jr P 



< 



1 



< 



< 



8 Jo 
1 



\\P 



\\p 



\\P 



7-3 



7-3 



7-3 



2 dr + C 



l dr + C 



2 dr 



5 




p 1 2 u y (j) 2 \u y y\p 1 2l dydr 



'0 JR 

.3 



<5 3 p 1+27 Jo 




Uyp 1 2 (fi 2 dydT 



R 




Uyp 1 2 (p 2 dydr. 



JR 



It follows from ( 13 .3 1 j) that 
h = 



'0 ^R £ 

< Ce5- 1 



< 



^j(j) y ip\dydr 




\p 2 7 u y (j) y ijj\dydT + Ce5 1 / / \p 1 1 u y (f) y 'ip\dydT 



\\P 



jr 

7-3 , 



2r-2 




'0 ./R 

(p-^ + p 1 - 37 )!^! 2 ^ 2 ^ 

I f Wp^^yfdr+CeH^pr^ [ T [ puyi, 2 dydr 

? _ „ ./O JR 



JR 



< 



< 



7-3 



\\p—<P v VdT + 




JR 



pu y ip 2 dydr, if e C 1. 



Finally, one has 




\ g mdydr<t I \\p^ y fdr + C 




lo Jr P 

Recalling that (13. 7p . ( 13 . 1 3[) and (I3.29p . one can get 

\g\ < \Uyy\ + \pU y 1p\ + C|p 7 " 2 p ; 

7—1 

< \Uyy\ + \pUy1P\ + C\p~Uy(j)\ 

Thus the last term in (13.381) can be estimated by 

( Li^ dydT 

u 2 
a yy 



JrP 



1+7 



dydr. 



< C 



< C 



< C 



< C 




lo JrP 

^3 



i^dydr + C 

1+7 a 




p^ul^dydr + C 



JR 




p u dydr 



/; 1 +7 / 
P Jo 



5p 1+ ^ 



<5 2 p 1+7 



r + -y 2 dr + C 



5p"> 
e 

5pT 




JR 




JR 



r JR 

{j)u y i\} 2 + p^^UyCp^dydr 
{jmyi\) 2 + Uyp 1 ^ '(ft 2 ) dydr 



+ - ( j (puyip 2 + Uyp 1 2 (p 2 )dydT, if e <C 1. 

° ^R 
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Substituting f l3T2?j) -f l3^0]) into f l3T26]) and recalling f[3TT2j) . it holds that 

+ p7~2 2 + 2 j dy+ ^ I ^2 + p7~2 -^2 + --^2 + -,-3^ ^ 



R v ^ «/ R 

< C7e (l/2-a)| lne |-l/2 ) if e< l. 

Step 3. As the last step, we estimate sup 

tG[0,ti] 



(3.41) 

y \\. For this, multiplying (I3.5P by 



-Tpyy/P S ives 



;f ), - + 4'» + «»f - "'(")^ + ^ = »t- (3 ' 42) 



Integrating the above equation over R 1 x [0, r] yields 



R z 



First, one has 




o JR 




'0 JR 



^yy<Py . *Pyy "0. 

P P 



jO 2 J 

(3.43) 




*/R 



p'(p)^ldydr 
P 



< 



\ff ^dydr + C Hip 
8 Jo Jr P Jo 



(()y\\ 2 dT. 



Then it follows from (EP0|) and that 



/o 7r P 




< 



1 

< - 




in 



yv dydr + C 




./R P 



-dydr 



lo JR P 

^dydr + Ce^-^e^ 1 ' 2 . 

Jr P 




Furthermore, we can compute that 



—dydr 



o ^R 




<C / ||^ w ||»||^||'dT 



< 



< - 



1 

< - 








:</;</ 



dydr + C 



3 dr 



f / —^-dydr + C sup HV'j/ll 1 

Jr P re[0,ri] 

-™dydr + C 

Jr P 



2 dr 




dr, 



where in the last inequality we have used the a priori assumption 



(3.44) 



(3.45) 



(3.46) 



SUP \\lpy\ 

re[o,n] 



< 1. 



(3.47) 
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Substituting ( j3HP , f l3^5|) and f l3^6|) into fl3T43]) and using f l37T3|) and (EHIjl . it holds that 
/ ^dy+ f / (n^ + ^^r^Ce^^llner 1 / 2 . (3.48) 

Therefore, (I3.17p can be derived directly from (I3.4ip and f!3.48j) and the a priori assumption 
(I3.47P is verified if e is suitably small. It follows from ( 13.1 7ft that if 1 < 7 < 2, then 

sup ||0(-,r)|| LO o <V2sup WH-.rW^Uyi;^ 2 

0<T<T! 0<T<n 

<C sup ( / (P-Vdy)H [ <P 2 y dyY (3-49) 

0<T<n JR Jr 

CCe^llnel" 1 / 4 , 



and if 7 > 2, then 



SUp \\<p{- : T L o 



0<T<Tl 0<T<T1 



< v/2 sup H^^r)!! 1 / 2 ^^-^)!! 1 / 2 



<C sup ( / fjP-*<fdy)*( I <p 2 y d y y 

0<r<n Jr Jr 



(3.50) 



< C/J ?e 4 2(7+4) |lne| 4 

1 1-7 

< Cei+ 4 \ hie 4 . 



And we also have 



sup \\i) y \\ < Ce { 3-f)|lne|- 1/4 < 1. (3.51) 

t6[0,ti] 

So from fl3.49p - fl3.5ip the a priori assumption (I3.10p is verified if e 1. On the other 
hand, by using Sobolev inequality, one can get 

sup WMWu- <V2 sup ^(•,r)|| 1 / 2 ||^(-,r)|| 1 /2 

0<T<T1 0<T<T! 

<V2^ sup IIVM-^iril^MU 1 / 2 

0<T<T! (^O.OZJ 

Ce 1/8 |lne|- 1/2 , if 1< 7 < 2, 



< 



<7 e 4(7+4) I lne|" 1/2 , if 7 > 2 



Thus the convergence rate (13.181) is justified and the proof of Lemma I3T21 is completed. 

□ 

Proof of Theorem ll.lt It remains to prove (11. 7p with a, b given in (I1.8P and 
respectively. From Lemma 12. 2\ Lemma 12.31 (iii) and Theorem 13.11 and recalling that 
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/i = e a \ lne|, 5 = e a , it holds that for any given positive constant h, there exist a constant 
Ch > which is independent of e such that 

sup||p(-,t) - p r (-)\\ L ™ 

t>h t 

sup ||0(-,t)|| £oo +sup||p(-,t) -^(i)IUoc +sup||p^(-) -ff{-)\\L°° 

I t>h. I I 



and 



< 



0<t<+oo t>h L t>h 

C h Ye 1/6 |lne|- 1/4 + 5|ln5| if 1< 7 < 2, 
C h ,(e^4|lne| (1 - 7)/4 + 5|ln5| + A if 7 > 2, 



< C h e a \lne\ 



sup ||m(-,t) - m r {-)\\ L oo 

t>h t 

< sup (\\m(-,t) -m(-,f)|| £ «» + ||m(-,t) -mj;(-)|| L « + ||m^(-) - m r (-) || £ „ 



< C sup (| 

0<T<+OO V 



+ sup ||m(-,t) - m r (-)\\ L ^ + \\m r (-) - m r (-)|| L o 
t>h y t t t 



< 



< 



CJe 1 / 8 ! lnel" 1 / 2 + e 1 / 6 ) hie)" 1 / 4 + 5| ln5| + fij , if 1< 7 < 2, 

C^eOT)|lne|- 1/2 + e^4|lne| (1 - 7)/4 + 5|ln5|+/i), if 7 > 2, 



C^llnel-a, if 1< 7 < 3, 
Cfte^l lne|, if 7 > 3, 



Thus the proof of Theorem 1.1 is completed. □ 
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